Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release.

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Computational Insight into Substrate-Induced Conformational Changes in Methionyl-tRNA Synthetase of Mycobacterium Tuberculosis.

Tuberculosis caused by Mycobacterium tuberculosis (M.tb) has killed millions worldwide. Antibiotic resistance leads to the ineffectiveness of the current therapies. Aminoacyl tRNA synthetase (aaRS) class of proteins involved in protein synthesis are promising bacterial targets for developing new therapies. Here, we carried out a systematic comparative study on the aaRS sequences from M.tb and human. We listed important M.tb aaRS that could be explored as potential M.tb targets alongside the detailed conformational space analysis of methionyl-tRNA synthetase (MetRS) in apo- and substrate-bound form, which is among the proposed targets. Understanding the conformational dynamics is central to the mechanistic understanding of MetRS, as the substrate binding leads to the conformational changes causing the reaction to proceed. We performed the most complete simulation study of M.tb MetRS for 6 microseconds (2 systems × 3 runs × 1 microsecond) in the apo and substrate-bound states. Interestingly, we observed differential features, showing comparatively large dynamics for the holo simulations, whereas the apo structures became slightly compact with reduced solvent exposed area. In contrast, the ligand size decreased significantly in holo structures possibly to relax ligand conformation. Our findings correlate with experimental studies, thus validating our protocol. Adenosine monophosphate moiety of the substrate exhibited quite higher fluctuations than the methionine. His21 and Lys54 were found to be the important residues forming prominent hydrogen bond and salt-bridge interactions with the ligand. The ligand-protein affinity decreased during simulations as computed by MMGBSA analysis over the last 500 ns trajectories, which indicates the conformational changes upon ligand binding. These differential features could be further explored for designing new M.tb inhibitors.

Identification of compounds contributing to umami taste of pea protein isolate.

The umami taste of pea protein ingredients can be desirable or undesirable based on the food application. The compounds contributing to the umami perception of pea protein isolate (PPI) were investigated. Sensory-guided prep-liquid chromatography fractionation of a 10% aqueous PPI solution revealed one well-known compound, monosodium glutamate (MSG), however, it was reported at a subthreshold concentration. Two umami enhancing compounds 5'-adenosine monophosphate (AMP) and 5'-uridine monophosphate (UMP) were subsequently identified after the LC fractions were re-evaluated with MSG. Sensory recombination studies, utilizing the aqueous PPI solution as the base, confirmed AMP and UMP were umami enhancers of MSG and contributed approximately 81% of the perceived umami intensity. However UMP was only reported to enhance umami perception in combination with AMP (not individually) indicating synergistic interactions were observed between the two enhancer compounds. Therefore the presence of all three compounds are important for umami perception and provide an improved basis to tailor the flavor profile in PPI products.

Energetic vs. entropic stabilization between a Remdesivir analogue and cognate ATP upon binding and insertion into the active site of SARS-CoV-2 RNA dependent RNA polymerase.

SARS-CoV-2 RNA dependent RNA polymerase (RdRp) serves as a highly promising antiviral drug target such as for a Remdesivir nucleotide analogue (RDV-TP or RTP). In this work, we mainly used alchemical all-atom simulations to characterize relative binding free energetics between the nucleotide analogue RTP and natural cognate substrate ATP upon initial binding and pre-catalytic insertion into the active site of SARS-CoV-2 RdRp. Natural non-cognate substrate dATP and mismatched GTP were also examined for computation control. We first identified significant differences in dynamical responses between nucleotide initial binding and subsequent insertion configurations to the open and closed active sites of the RdRp, respectively, though the RdRp protein conformational changes between the active site's open and closed states are subtle. Our alchemical simulations indicated that upon initial binding (active site open), RTP and ATP show similar binding free energies to the active sites while in the insertion state (active site closed), ATP is more stabilized (∼-2.4 kcal mol-1) than RTP in free energetics. Additional analyses show, however, that RTP is more stabilized in binding energetics than ATP, in both the insertion and initial binding states, with RTP more stabilized due to the electrostatic energy in the insertion state and due to vdW energy in the initial binding state. Hence, it appears that natural cognate ATP still excels at association stability with the RdRp active site due to that ATP maintains sufficient flexibilities e.g., in base pairing with the template, which exemplifies an entropic contribution to the cognate substrate stabilization. These findings highlight the importance of substrate flexibilities in addition to energetic stabilization in antiviral nucleotide analogue design.

Guvermectin Biosynthesis Revealing the Key Role of a Phosphoribohydrolase and Structural Insight into the Active Glutamate of a Noncanonical Adenine Phosphoribosyltransferase.

Guvermectin is a novel plant growth regulator that has been registered as a new agrochemical in China. It is an adenosine analogue with an unusual psicofuranose instead of ribose. Herein, the gene cluster responsible for guvermectin biosynthesis in Streptomyces caniferus NEAU6 is identified using gene interruption and heterologous expression experiments. A key intermediate psicofuranine 6'-phosphate (PMP) is chemically synthesized, and the functions of GvmB, C, D, and E are verified by individual stepwise enzyme reactions in vitro. The results also show that the biosynthesis of guvermectin is coupled with adenosine production by a single cluster. The higher catalytic efficiency of GvmB on PMP than AMP ensures the effective biosynthesis of guvermectin. Moreover, a phosphoribohydrolase GvmA is employed in the pathway that can hydrolyze AMP but not PMP and shows higher catalytic efficiency for the AMP hydrolysis than that of the AMP dephosphorylation by GvmB, leading to shunting of adenosine biosynthesis toward the production of guvermectin. Finally, the crystal structure of GvmE in complex with the product PMP has been solved. Glu160 at the C-terminal is identified as the acid/base for protonation/deprotonation of N7 of the adenine ring, demonstrating that GvmE is a noncanonical adenine phosphoribosyltransferase.

The binding model of adenosine-specific DNA aptamer: Umbrella sampling study.

We report a novel model of the selective binding mechanism of adenosine-specific DNA aptamer. Our theoretical investigations of AMP (Adenosine monophosphate) dissociation from aptamer-AMP complexes reveals new details of aptamer molecular specificity and stabilisation factors. Umbrella sampling MD calculations using parmbsc1 force field shows that the disordered structure of the internal loop of the unbound aptamer hairpin has a characteristic packing of guanines, which prevents barrier-free penetration of ligands into the site cavity. Also, this disordered structure of the unbound aptamer has a network of hydrogen bonds stabilising the cavity near the target guanines within the binding sites during the whole binding process. We suggested that the first AMP molecule binds to the disordered structure of the site closest to the aptamer hairpin stem and spends some free energy on ordering of the internal loop. Then the second AMP molecule binds to the ordered site closest to the aptamer hairpin loop with a lower energy gain. As a result, the induced-fit binding model is the most applicable for this aptamer and does not contradict the modern experimental NMR and calorimetry data.

Insights into the binding of arginine to adenosine phosphate from mimetic complexes.

The amino acid arginine plays a key role in the interaction of proteins with adenosine phosphates, as its protonated guanidinium side group is capable of building multipodal H-bonding interactions with the oxygen atoms of the phosphate, phosphoester and ribose moieties and with the nitrogen atoms of adenine. Protein interactions often take place in competition with other ionic species, typically metal cations, which are prone to build concerted coordination arrangements with the same centers of negative charge as guanidinium. We report on a vibrational spectroscopy and computational investigation of a positively charged ternary complex formed by adenosine monophosphate (AMP) with methyl guanidinium and Na+. Following a bottom-up approach, an analogous complex with ribose phosphate is characterized as well, which serves to assess the individual role of the phosphate, sugar and adenine moieties in the binding process and to compare, within a single complex, the interactions associated with diffuse versus localized charge distributions of guanidinium and the alkali cation, respectively. The results indicate that Na+ is preferentially hosted in a semi-rigid pocket formed by the phosphoester-adenosine backbone of AMP and displaces guanidinium to a peripheral binding to the phosphate anionic end group. This suggests that the control of the salt concentration may constitute an effective route to modulate protein-AMP complexation.

Direct adenylation from 5'-OH-terminated oligonucleotides by a fusion enzyme containing Pfu RNA ligase and T4 polynucleotide kinase.

5'-Adenylated oligonucleotides (AppOligos) are widely used for single-stranded DNA/RNA ligation in next-generation sequencing (NGS) applications such as microRNA (miRNA) profiling. The ligation between an AppOligo adapter and target molecules (such as miRNA) no longer requires ATP, thereby minimizing potential self-ligations and simplifying library preparation procedures. AppOligos can be produced by chemical synthesis or enzymatic modification. However, adenylation via chemical synthesis is inefficient and expensive, while enzymatic modification requires pre-phosphorylated substrate and additional purification. Here we cloned and characterized the Pfu RNA ligase encoded by the PF0353 gene in the hyperthermophilic archaea Pyrococcus furiosus. We further engineered fusion enzymes containing both Pfu RNA ligase and T4 polynucleotide kinase. One fusion enzyme, 8H-AP, was thermostable and can directly catalyze 5'-OH-terminated DNA substrates to adenylated products. The newly discovered Pfu RNA ligase and the engineered fusion enzyme may be useful tools for applications using AppOligos.

In Vitro Selection of Remdesivir-Resistant SARS-CoV-2 Demonstrates High Barrier to Resistance.

In vitro selection of remdesivir-resistant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) revealed the emergence of a V166L substitution, located outside of the polymerase active site of the Nsp12 protein, after 9 passages of a single lineage. V166L remained the only Nsp12 substitution after 17 passages (10 µM remdesivir), conferring a 2.3-fold increase in 50% effective concentration (EC50). When V166L was introduced into a recombinant SARS-CoV-2 virus, a 1.5-fold increase in EC50 was observed, indicating a high in vitro barrier to remdesivir resistance.

Drug-Membrane Interactions: Effects of Virus-Specific RNA-Dependent RNA Polymerase Inhibitors Remdesivir and Favipiravir on the Structure of Lipid Bilayers.

The two RNA-dependent RNA polymerase inhibitors remdesivir and favipiravir were originally developed and approved as broad-spectrum antiviral drugs for the treatment of harmful viral infections such as Ebola and influenza. With the outbreak of the global SARS-CoV-2 pandemic, the two drugs were repurposed for the treatment of COVID-19 patients. Clinical studies suggested that the efficacy of the drugs is enhanced in the case of an early or even prophylactic application. Because the contact between drug molecules and the plasma membrane is essential for a successful permeation process of the substances and therefore for their intracellular efficiency, drug-induced effects on the membrane structure are likely and have already been shown for other substances. We investigated the impact of remdesivir and favipiravir on lipid bilayers in model and cell membranes via several biophysical approaches. The measurements revealed that the embedding of remdesivir molecules in the lipid bilayer results in a disturbance of the membrane structure of the tested phospholipid vesicles. Nevertheless, in a cell-based assay, the presence of remdesivir induced only weak hemolysis of the treated erythrocytes. In contrast, no experimental indication for an effect on the structure and integrity of the membrane was detected in the case of favipiravir. Regarding potential prophylactic or accompanying use of the drugs in the therapy of COVID-19, the physiologically relevant impacts associated with the drug-induced structural modifications of the membrane might be important to understand side effects and/or low effectivities.

ADP-Induced Conformational Transition of Human Adenylate Kinase 1 Is Triggered by Suppressing Internal Motion of α3α4 and α7α8 Fragments on the ps-ns Timescale.

Human adenylate kinase 1 (hAK1) plays a vital role in the energetic and metabolic regulation of cell life, and impaired functions of hAK1 are closely associated with many diseases. In the presence of Mg2+ ions, hAK1 in vivo can catalyze two ADP molecules into one ATP and one AMP molecule, activating the downstream AMP signaling. The ADP-binding also initiates AK1 transition from an open conformation to a closed conformation. However, how substrate binding triggers the conformational transition of hAK1 is still unclear, and the underlying molecular mechanisms remain elusive. Herein, we determined the solution structure of apo-hAK1 and its key residues for catalyzing ADP, and characterized backbone dynamics characteristics of apo-hAK1 and hAK1-Mg2+-ADP complex (holo-hAK1) using NMR relaxation experiments. We found that ADP was primarily bound to a cavity surrounded by the LID, NMP, and CORE domains of hAK1, and identified several critical residues for hAK1 catalyzing ADP including G16, G18, G20, G22, T39, G40, R44, V67, D93, G94, D140, and D141. Furthermore, we found that apo-hAK1 adopts an open conformation with significant ps-ns internal mobility, and Mg2+-ADP binding triggered conformational transition of hAK1 by suppressing the ps-ns internal motions of α3α4 in the NMP domain and α7α8 in the LID domain. Both α3α4 and α7α8 fragments became more rigid so as to fix the substrate, while the catalyzing center of hAK1 experiences promoted µs-ms conformational exchange, potentially facilitating catalysis reaction and conformational transition. Our results provide the structural basis of hAK1 catalyzing ADP into ATP and AMP, and disclose the driving force that triggers the conformational transition of hAK1, which will deepen understanding of the molecular mechanisms of hAK1 functions.

A Similarity-Based Method for Predicting Enzymatic Functions in Yeast Uncovers a New AMP Hydrolase.

Despite decades of research and the availability of the full genomic sequence of the baker's yeast Saccharomyces cerevisiae, still a large fraction of its genome is not functionally annotated. This hinders our ability to fully understand cellular activity and suggests that many additional processes await discovery. The recent years have shown an explosion of high-quality genomic and structural data from multiple organisms, ranging from bacteria to mammals. New computational methods now allow us to integrate these data and extract meaningful insights into the functional identity of uncharacterized proteins in yeast. Here, we created a database of sensitive sequence similarity predictions for all yeast proteins. We use this information to identify candidate enzymes for known biochemical reactions whose enzymes are unidentified, and show how this provides a powerful basis for experimental validation. Using one pathway as a test case we pair a new function for the previously uncharacterized enzyme Yhr202w, as an extra-cellular AMP hydrolase in the NAD degradation pathway. Yhr202w, which we now term Smn1 for Scavenger MonoNucleotidase 1, is a highly conserved protein that is similar to the human protein E5NT/CD73, which is associated with multiple cancers. Hence, our new methodology provides a paradigm, that can be adopted to other organisms, for uncovering new enzymatic functions of uncharacterized proteins.

Synthesis of Nucleoside-like Molecules from a Pyrolysis Product of Cellulose and Their Computational Prediction as Potential SARS-CoV-2 RNA-Dependent RNA Polymerase Inhibitors.

(1R,5S)-1-Hydroxy-3,6-dioxa-bicyclo[3.2.1]octan-2-one, available by an efficient catalytic pyrolysis of cellulose, has been applied as a chiral building block in the synthesis of seven new nucleoside analogues, with structural modifications on the nucleobase moiety and on the carboxyl- derived unit. The inverted configuration by Mitsunobu reaction used in their synthesis was verified by 2D-NOESY correlations, supported by the optimized structure employing the DFT methods. An in silico screening of these compounds as inhibitors of SARS-CoV-2 RNA-dependent RNA polymerase has been carried out in comparison with both remdesivir, a mono-phosphoramidate prodrug recently approved for COVID-19 treatment, and its ribonucleoside metabolite GS-441524. Drug-likeness prediction and data by docking calculation indicated compound 6 [=(3S,5S)-methyl 5-(hydroxymethyl)-3-(6-(4-methylpiperazin-1-yl)-9H-purin-9-yl)tetrahydrofuran-3-carboxylate] as the best candidate. Furthermore, molecular dynamics simulation showed a stable interaction of structure 6 in RNA-dependent RNA polymerase (RdRp) complex and a lower average atomic fluctuation than GS-441524, suggesting a well accommodation in the RdRp binding pocket.

Potential COVID-19 Drug Candidates Based on Diazinyl-Thiazol-Imine Moieties: Synthesis and Greener Pastures Biological Study.

A novel series of 1-aryl-N-[4-phenyl-5-(arylazo)thiazol-2-yl)methanimines has been synthesized via the condensation of 2-amino-4-phenyl-5-arylazothiazole with various aromatic aldehydes. The synthesized imines were characterized by spectroscopic techniques, namely 1H and 13C-NMR, FTIR, MS, and Elemental Analysis. A molecular comparative docking study for 3a-f was calculated, with reference to two approved drugs, Molnupiravir and Remdesivir, using 7BQY (Mpro; PDB code 7BQY; resolution: 1.7 A°) under identical conditions. The binding scores against 7BQY were in the range of -7.7 to -8.7 kcal/mol for 3a-f. The high scores of the compounds indicated an enhanced binding affinity of the molecules to the receptor. This is due to the hydrophobic interactions and multi-hydrogen bonds between 3a-f ligands and the receptor's active amino acid residues. The main aim of using in silco molecular docking was to rank 3a-f with respect to the approved drugs, Molnupiravir and Remdesivir, using free energy methods as greener pastures. A further interesting comparison presented the laydown of the ligands before and after molecular docking. These results and other supporting statistical analyses suggested that ligands 3a-f deserve further investigation in the context of potential therapeutic agents for COVID-19. Free-cost, PASS, SwissADME, and Way2drug were used in this research paper to determine the possible biological activities and cytotoxicity of 3a-f.

Efficient incorporation and template-dependent polymerase inhibition are major determinants for the broad-spectrum antiviral activity of remdesivir.

Remdesivir (RDV) is a direct-acting antiviral agent that is approved in several countries for the treatment of coronavirus disease 2019 caused by the severe acute respiratory syndrome coronavirus 2. RDV exhibits broad-spectrum antiviral activity against positive-sense RNA viruses, for example, severe acute respiratory syndrome coronavirus and hepatitis C virus, and nonsegmented negative-sense RNA viruses, for example, Nipah virus, whereas segmented negative-sense RNA viruses such as influenza virus or Crimean-Congo hemorrhagic fever virus are not sensitive to the drug. The reasons for this apparent efficacy pattern are unknown. Here, we expressed and purified representative RNA-dependent RNA polymerases and studied three biochemical parameters that have been associated with the inhibitory effects of RDV-triphosphate (TP): (i) selective incorporation of the nucleotide substrate RDV-TP, (ii) the effect of the incorporated RDV-monophosphate (MP) on primer extension, and (iii) the effect of RDV-MP in the template during incorporation of the complementary UTP. We found a strong correlation between antiviral effects and efficient incorporation of RDV-TP. Inhibition in primer extension reactions was heterogeneous and usually inefficient at higher NTP concentrations. In contrast, template-dependent inhibition of UTP incorporation opposite the embedded RDV-MP was seen with all polymerases. Molecular modeling suggests a steric conflict between the 1'-cyano group of the inhibitor and residues of the structurally conserved RNA-dependent RNA polymerase motif F. We conclude that future efforts in the development of nucleotide analogs with a broader spectrum of antiviral activities should focus on improving rates of incorporation while capitalizing on the inhibitory effects of a bulky 1'-modification.

Genomic diversity and molecular dynamics interaction on mutational variances among RB domains of SARS-CoV-2 interplay drug inactivation.

The scientific community has been releasing whole genomic sequences of SARS-CoV-2 to facilitate the investigation of molecular features and evolutionary history. We retrieved 36 genomes of 18 prevalent countries of Asia, Europe and America for genomic diversity and mutational analysis. Besides, we studied mutations in the RBD regions of Spike (S) proteins to analyze the drug efficiency against these mutations. In this research, phylogenenetic analysis, evolutionary modeling, substitution pattern analysis, molecular docking, dynamics simulation, etc. were performed. The genomic sequences showed >99% similarity with the reference sequence of China.TN93 + G was predicted as a best nucleotide substitution model. It was revealed that effective transition from the co-existing SARS genome to the SARS-CoV-2 and a noticeable positive selection in the SARS-CoV-2 genomes occurred. Moreover, three mutations in RBD domain, Val/ Phe367, Val/ Leu 382 and Ala/ Val522, were discovered in the genomes from Netherland, Bangladesh and the USA, respectively. Molecular docking and dynamics study showed RBD with mutation Val/Leu382 had the lowest binding affinity with remdesivir. In conclusion, the SARS-CoV-2 genomes are similar, but multiple degrees of transitions and transversions occurred. The mutations cause a significant conformational change, which are needed to be investigated during drug and vaccine development.

Longitudinal analysis of SARS-CoV-2 spike and RNA-dependent RNA polymerase protein sequences reveals the emergence and geographic distribution of diverse mutations.

Amid the ongoing COVID-19 pandemic, it has become increasingly important to monitor the mutations that arise in the SARS-CoV-2 virus, to prepare public health strategies and guide the further development of vaccines and therapeutics. The spike (S) protein and the proteins comprising the RNA-Dependent RNA Polymerase (RdRP) are key vaccine and drug targets, respectively, making mutation surveillance of these proteins of great importance. Full protein sequences were downloaded from the GISAID database, aligned, and the variants identified. 437,006 unique viral genomes were analyzed. Polymorphisms in the protein sequence were investigated and examined longitudinally to identify sequence and strain variants appearing between January 5th, 2020 and January 16th, 2021. A structural analysis was also performed to investigate mutations in the receptor binding domain and the N-terminal domain of the spike protein. Within the spike protein, there were 766 unique mutations observed in the N-terminal domain and 360 in the receptor binding domain. Four residues that directly contact ACE2 were mutated in more than 100 sequences, including positions K417, Y453, S494, and N501. Within the furin cleavage site of the spike protein, a high degree of conservation was observed, but the P681H mutation was observed in 10.47% of sequences analyzed. Within the RNA dependent RNA polymerase complex proteins, 327 unique mutations were observed in Nsp8, 166 unique mutations were observed in Nsp7, and 1157 unique mutations were observed in Nsp12. Only 4 sequences analyzed contained mutations in the 9 residues that directly interact with the therapeutic Remdesivir, suggesting limited mutations in drug interacting residues. The identification of new variants emphasizes the need for further study on the effects of the mutations and the implications of increased prevalence, particularly for vaccine or therapeutic efficacy.

Direct and Catalytic C-Glycosylation of Arenes: Expeditious Synthesis of the Remdesivir Nucleoside.

Since early 2020, scientists have strived to find an effective solution to fight SARS-CoV-2, in particular by developing reliable vaccines that inhibit the spread of the disease and repurposing drugs for combatting its effects on the human body. The antiviral prodrug Remdesivir is still the most widely used therapeutic during the early stages of the infection. However, the current synthetic routes rely on the use of protecting groups, air-sensitive reagents, and cryogenic conditions, thus impeding a cost-efficient supply to patients. We have, therefore, focused on the development of a straightforward, direct addition of (hetero)arenes to unprotected sugars. Here we report a silylium-catalyzed and completely stereoselective C-glycosylation that initially yields the open-chain polyols, which can be selectively cyclized to provide either the kinetic α-furanose or the thermodynamically favored ß-anomer. The method significantly expedites the synthesis of Remdesivir precursor GS-441524 after a subsequent Mn-catalyzed C-H oxidation and deoxycyanation.

Human carboxylesterase 1A plays a predominant role in the hydrolytic activation of remdesivir in humans.

Remdesivir, an intravenous nucleotide prodrug, has been approved for treating COVID-19 in hospitalized adults and pediatric patients. Upon administration, remdesivir can be readily hydrolyzed to form its active form GS-441524, while the cleavage of the carboxylic ester into GS-704277 is the first step for remdesivir activation. This study aims to assign the key enzymes responsible for remdesivir hydrolysis in humans, as well as to investigate the kinetics of remdesivir hydrolysis in various enzyme sources. The results showed that remdesivir could be hydrolyzed to form GS-704277 in human plasma and the microsomes from human liver (HLMs), lung (HLuMs) and kidney (HKMs), while the hydrolytic rate of remdesivir in HLMs was the fastest. Chemical inhibition and reaction phenotyping assays suggested that human carboxylesterase 1 (hCES1A) played a predominant role in remdesivir hydrolysis, while cathepsin A (CTSA), acetylcholinesterase (AchE) and butyrylcholinesterase (BchE) contributed to a lesser extent. Enzymatic kinetic analyses demonstrated that remdesivir hydrolysis in hCES1A (SHUTCM) and HLMs showed similar kinetic plots and much closed Km values to each other. Meanwhile, GS-704277 formation rates were strongly correlated with the CES1A activities in HLM samples from different individual donors. Further investigation revealed that simvastatin (a therapeutic agent for adjuvant treating COVID-19) strongly inhibited remdesivir hydrolysis in both recombinant hCES1A and HLMs. Collectively, our findings reveal that hCES1A plays a predominant role in remdesivir hydrolysis in humans, which are very helpful for predicting inter-individual variability in response to remdesivir and for guiding the rational use of this anti-COVID-19 agent in clinical settings.

Structural Bioinformatics Used to Predict the Protein Targets of Remdesivir and Flavones in SARS-CoV-2 Infection.

BACKGROUND: During the current SARS-CoV-2 pandemic, the identification of effective antiviral drugs is crucial. Unfortunately, no specific treatment or vaccine is available to date. OBJECTIVE: Here, we aimed to predict the interactions with SARS-CoV-2 proteins and protein targets from the human body for some flavone molecules (kaempferol, morin, pectolinarin, myricitrin, and herbacetin) in comparison to synthetic compounds (hydroxychloroquine, remdesivir, ribavirin, ritonavir, AMD-070, favipiravir). METHODS: Using MOE software and advanced bioinformatics and cheminformatics portals, we conducted an extensive analysis based on various structural and functional features of compounds, such as their amphiphilic field, flexibility, and steric features. The structural similarity analysis of natural and synthetic compounds was performed using Tanimoto coefficients. The interactions of some compounds with SARS-CoV-2 3CLprotease or RNA-dependent RNA polymerase were described using 2D protein-ligand interaction diagrams based on known crystal structures. The potential targets of considered compounds were identified using the SwissTargetPrediction web tool. RESULTS: Our results showed that remdesivir, pectolinarin, and ritonavir present a strong structural similarity which may be correlated to their similar biological activity. As common molecular targets of compounds in the human body, ritonavir, kaempferol, morin, and herbacetin can activate multidrug resistance-associated proteins, while remdesivir, ribavirin, and pectolinarin appear as ligands for adenosine receptors. CONCLUSION: Our evaluation recommends remdesivir, pectolinarin, and ritonavir as promising anti- SARS-CoV-2 agents.